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We examine whether the annual modulation found by the DAMA dark matter experiment can be 
explained by Weakly Interacting Massive Particles ( WIMPs) , in light of new null results from other 
experiments. CDMS II has already ruled out most WIMP-nucleus spin-independent couplings as an 
explanation for DAMA data. Hence we here focus on spin-dependent (axial vector; SD) couplings 
of WIMPs to nuclei. We expand upon previous work by (i) considering the general case of coupling 
to both protons and neutrons and (ii) incorporating bounds from all existing experiments. We 
note the surprising fact that CMDS II places one of the strongest bounds on the WIMP-neutron 
cross-section, and show that SD WIMP-neutron scattering alone is excluded. We also show that 
SD WIMP-proton scattering alone is allowed only for WIMP masses in the 5-13 GeV range. For 
the general case of coupling to both protons and neutrons, we find that, for WIMP masses above 
13 GeV and below 5 GeV, there is no region of parameter space that is compatible with DAMA 
and all other experiments. In the range (5-13) GeV, we find acceptable regions of parameter space, 
including ones in which the WIMP-neutron coupling is comparable to the WIMP-proton coupling. 



I. INTRODUCTION 

The dark halo of our Galaxy may consist of WIMPs 
(Weakly Interacting Massive Particles). Numerous col- 
laborations worldwide have been searching for these par- 
ticles. Direct detection experiments attempt to observe 
the nuclear recoil caused by these dark matter particles 
interacting with nuclei in the detectors. Indirect detec- 
tion experiments search for the signatures of WIMP an- 
nihilation in the Sun, the Earth, or the Galactic Halo. It 
is well known that the count rate in WIMP direct detec- 
tion experiments will experience an annual modulation 
P, Q as a result of the motion of the Earth around the 
Sun: the relative velocity of the detector with respect to 
the WIMPs depends on the time of year. 

The discovery of an annual modulation by the 
DAMA/Nal experiment (hereafter, "DAMA") is the 
only positive signal seen in any dark matter search. How- 
ever, recent null results from other experiments, particu- 
larly EDELWEISS and CDMS II, severely bound the pa- 
rameter space of WIMPs that could possibly explain the 
DAMA data. In fact, if the interactions between WIMPs 
and nuclei are spin- independent (SI), then CDMS II has 
ruled out DAMA (with some exceptions 0). In this pa- 
per, we will consider the alternative of spin-dependent 
(SD) WIMP cross-sections. We expand upon previous 
work by (i) considering the general case of coupling to 
both protons and neutrons and (ii) incorporating bounds 
from all existing experiments, including, in particular, 
the null results recently released by CDMS II and Super- 
Kamiokandc (Super-K), which were not available to pre- 
vious authors. Previously, UUio et al. studied the cases 
of interactions with protons only or neutrons only Q 
and made a comparison with experiments that had been 
performed at that time; they argued that some cases of 
mixed interactions would be ruled out as well. In 2003 



Kurylov and Kamionkowski performed a general analy- 
sis for point-like WIMPs of arbitrary spin to find what 
WIMPs were compatible with DAMA and the null exper- 
iments at the time Giuliani looked at the more gen- 
eral case of mixed interactions , but only applied the 
analysis to DAMA for a limited choice of WIMP mass. 
We here perform a general analysis with coupling to 
both protons and neutrons to see if any parameter space 
remains that could explain simultaneously the DAMA 
and CDMS data. To be complete, we will also incorpo- 
rate bounds from other experiments, including Super-K 
(which measures indirect detection in the Sun). Notice, 
however, that these indirect detection limits do not ap- 
ply if WIMPs do not effectively annihilate in the Sun, 
as for example if the WIMP is not its own anti-particle 
and there is a sufficiently strong cosmic asymmetry in 
the number of WIMPs and anti- WIMPs. 

We find those regions of WIMP parameter space that 
survive all existing experiments for the case of spin- 
dependent interactions. Although the current CDMS 
data are more stringent than those studied in previous 
work, our final bounds are somewhat less restrictive for 
WIMP-proton coupling than those quoted by UUio et al. 
because the latest results released by Super-K are less 
restrictive than those approximated by UUio et al. and 
because we are reluctant to extend Super-K data beyond 
those published by the experimentalists (as elaborated in 
the Discussion section at the end of the paper). 

We also note that, although CDMS data have typi- 
cally been ignored in the spin-dependent sector due to 
the small natural abundance of spin odd isotopes in the 
detector, in fact CDMS II places one of the strongest 
bounds on the WIMP-neutron cross-section {n.h. even 
CDMS I placed interesting bounds particularly at low 
masses) . 
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II. 



WIMP DETECTION 



The basic idea underlying WIMP direct detection is 
straightforward: the experiments seek to measure the en- 
ergy deposited when a WIMP interacts with a nucleus in 
the detector 

If a WIMP of mass m scatters elastically from a nu- 
cleus of mass M, it will deposit a recoil energy E = 
{li^v'^IM){l - cos 6*), where ^ = mM/{m + M) is the 
reduced mass, v is the speed of the WIMP relative to 
the nucleus, and 9 is the scattering angle in the center of 
mass frame. The differential recoil rate per unit detec- 
tor mass for a WIMP mass m, typically given in units of 
counts/kg/day/keV, can be written as: 
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where p = 0.3 GeV/cm^ is the standard local halo WIMP 
density, q — V 2ME is the nucleus recoil momentum, 
a{q) is the WIMP-nucleus cross-section, and information 
about the WIMP velocity distribution is encoded into the 
mean inverse speed rj{E,t), 
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represents the minimum WIMP velocity that can result 
in a recoil energy E and /d(u, is the (usually time- 
dependent) distribution of WIMP velocities u relative to 
the detector. 

For WIMPs in the Milky Way halo, the most frequently 
employed WIMP velocity distribution is that of a sim- 
ple isothermal sphere In such a model, the Galactic 
WIMP speeds with respect to the halo obey a Maxwellian 
distribution with a velocity dispersion ah truncated at 
the escape velocity Uosc, 



It is well known that the count rate in WIMP detectors 
will experience an annual modulation P, Q as a result 
of the motion of the Earth around the Sun: the rela- 
tive velocity of the detector with respect to the WIMPs 
depends on the time of year. The disk of the galaxy ro- 
tates through the relatively stationary halo, giving the 
Sun a velocity = 232 km/s relative to the halo. In 
addition, the Earth travels about the Sun at = 30 
km/s (relative to the Sun), with the Sun's motion in the 
Galaxy being at 60deg from the Earth's orbital plane 
(v© = v© V©). 

We will use the isothermal halo model in our calcula- 
tions, but its validity is by no means guaranteed - the 
actual halo may be squashed or anisotropic , may con- 
tain streams from smaller galaxies in the process of being 
absorbed by the Milky Way |11J or may, in fact, not yet 
be thermalized. 

An alternative way to search for WIMP dark matter of 
relevance to this paper is via indirect detection of WIMP 
annihilation in the Sun. When WIMPs pass throu gh a 
large celestial body, such as the Sun or the Earth [Hill, 
interactions can lead to gravitational capture if enough 
energy is lost in the collision to fall below the escape ve- 
locity. Captured WIMPs soon fall to the body's core and 
eventually annihilate with other WIMPs. These annihila- 
tions lead to high-energy neutrinos that can be observed 
by Earth-based detectors such as Supcr-Kamiokande [l4l . 
Baksan [l^, and AMANDA O, and the IceCube \l% 
and ANTARES (T3| projects. The annihilation rate de- 
pends on the capture rate of WIMPs, which is in turn de- 
termined by the WIMP scattering cross section off nuclei 
in the celestial body. While the Earth is predominantly 
composed of spinless nuclei, the Sun is mostly made of 
hydrogen, which has spin. Thus the spin-dependent cross 
section of WIMPs off nucleons can be probed by measur- 
ing the annihilation signals from WIMP annihilation in 
the Sun. Other indirect detection methods search for 
WIMPs that annihilate in the Galactic Halo or near the 
Galactic Center where they produce neutrinos, positrons, 
or antiprotons that may be seen in detectors on the Earth 
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Ncsc — erf(z) — 2zexp(— 2;^)/7r 
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with z — Wcsc/i^Oj is a normalization factor. For the sake 
of illustration, we take ah = 270 km/s and Vcsc — 650 
km/s. The results for the mean inverse speed ri{E,t) for 
the case of the isothermal sphere haye_previously been 
calculated and can be found e.g. in 



III. 



SPIN-INDEPENDENT CROSS-SECTION 



The cross-section for spin-independent WIMP interac- 
tions is given by 



(7) 



where ctq is the zero-momentum WIMP-nuclear cross- 
section and F{q) is the nuclear form factor, normalized 
to F{0) = 1 ; a description of these form factors may be 
found in [2J and [IJ. For purely scalar interactions, 



f^O, scalar 



[Zfp + {A-Z)U 



(8) 



Here Z is the number of protons, ^4 — Z is the number 
of neutrons, and fp and /„ are the WIMP couplings to 
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the proton and nucleon, respectively. In most instances, 



fp] the WIMP-nucleus cross-section can then be 



fn 

given in terms of the WIMP-proton cross-section as a 
result of Eqn. (|SJ): 



A' 



(9) 



where the fip is the proton- WIMP reduced mass, and A 
is the atomic mass of the target nucleus. 

In this model, for a given WIMP mass, tXp is the only 
free parameter and its limits are easily calculated for a 
given experiment. These limits are routinely calculated 
by the dark matter experiments and it has been found 
that the DAMA and CDMS results are incompatible for 
most WIMP mass in the spin-independent case, assuming 
the standard isothermal halo [23 (an exception is the 
case of light WIMPs of 6-9 GeV discussed in 



IV. 



SPIN-DEPENDENT CROSS-SECTION 



For the remainder of the paper, we will focus on 
spin-dependent interactions. The generic form for 
the spin-dependent WIMP-nucleus cross-section includes 
two couplings- the WIMP-proton coupling Op and the 
WIMP-neutron coupling a„. 
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Here the nuclear structure functions Spp{q), Snn{q), and 
Spn{q) are functions of the exchange momentum q and 
are specific to each nucleus. The quantities ap and a„ are 
actually the axial four-fermion WIMP-nucleon couplings 
in units of 2y/2GF [Ulllil- 

For a given experiment and WIMP mass, one can inte- 
grate Eqn. 1^ over an energy bin E1-E2, using the above 
cross-section and factoring in efficiencies, quenching fac- 
tors, etc., to obtain the expected number of recoil events 
Nrec as a function of the parameters ap and a„: 

Nrec(^0,p, (In) = ^reclp + Bj-ecClpCln + CrecO-n- (Ha) 
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rec; and Cfec 



are constants that can be calculated 



from the integration (see Appendix) . They differ between 
experiments and they depend on the WIMP mass and 
WIMP velocity distribution. The expected value of the 
amplitude of the annual modulation can be similarly cal- 
culated as Eqn. (|lla|l : 



A^ma(^p: ^n) — ^ma'^p ^ ^ma^p^n ~t- 



(lib) 



where N^a is the modulation amplitude and Ama, Bma, 
and Cma will again be calculable strictly from the ex- 
perimental parameters and the WIMP mass and velocity 
distribution. A more detailed discussion of calculating 
Nreci Nma, and the associated constants is given in the 



Appendix. Null search results place an upper limit on 
Nrec {Nrec < Nrec,max), which then constrains the pos- 



sible values for ap and a 



A positive modulation signal, 

mm ^ ^7na ^ ^ma,max ) : likcwisC 



such as in DAMA {N^a 
constrains the possible values for Op and a„. 

Previous authors, when searching for regions compat- 
ible with both the DAMA signal and the null results of 
other experiments , have studied the following two spe- 
cific cases: (i) a„ = so that only protons contribute to 
WIMP interactions, and (ii) ap = so that only neu- 
trons contribute. Here we examine the more general case 
where both ap and a„ may be non-zero (some of this pa- 
rameter space has also been examined by Giuliani j2j). 
To do this, we examine the form of the limits on ap and 
a„ implied by Eqn. (|11|) . 

Eqn. Hll|) describes a conic section in the ap-a„ plane. 
This conic section can be an ellipse, a hyperbola, or a 
set of two straight lines (it cannot be a parabola because 
linear terms in Op or a„ are absent). Which is the case 
depends on the values of the coefficients A, B, and C. 
For the benefit of the reader, we will briefly review the 
relevant shapes here (our conies are always centered at 
the origin, so we only describe such cases below). 

An ellipse in the x-y plane can be written as: 



62 



(12) 



This ellipse is symmetric about the x and y axes. Adding 
a cross-term (i.e. xy) in Eqn. (|12l) (like the B term in 
Eqn. (Illfl ) essentially rotates the ellipse- the major and 
minor axes do not lie on the x and y axes, but along some 
new rotated axes. 

Another conic section is the hyperbola: 
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= 1 



(13) 



which, in this form, gives two branches of a hyperbola 
open along the -fa; and —x directions. Rotation of this 
form also leads to a cross-term. 

Two parallel lines are also a conic (technically, a "de- 
generate" conic), that can be written in the form, e.g.: 



x' - a' (14) 

which describes two lines parallel to the y-axis, located at 
X = +a and x = —a. Other parallel lines can be rotated 
into this form as well. 

Note that, while other conies exist besides the three 
described, these are the only forms possible for Eqn. Hll|l 
(which is not the most general second degree polyno- 
mial). 

Eqn. Hll|l can, under suitable rotation in the ap-a„ 
plane, be expressed in one of the forms in Eqns. H12|l . 
(|13|l . or H14|) . The shape of our conic can be determined 
by the value of the discriminant 
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Experiment 


Exposure [kg-day] 


iiiresiioia [keVJ 


rjiticiency [/o\ 


Constraint (for stated recoil energies) 


Kei. 


OUlVlb i 


cc . a ceo 
Bl: D.08 





E <10 keV: 7.6 
£ < 20 keV: 22.8 
£ > 20 keV: 38 


0—00 keV: <2.o events [j) 


I28| 




^e. oz.D 


1 n 
iU 


iU— iUU Kev. <z.o events (^jj 




CDMS ii 
(projection) 


Si: 300 
Ge: 1200 


5 


5-100 keV: <2.3 events (f) 


[29] 


iltDtliLiW tlilOu 




on 


iUU 


on inn ^o tiiinan+c /'+\ 
zu~iuu Kev. <.z.o events (^jj 


loU] 


CRESST I 


AI2O3: 1.51 


0.6 


100 


(t) 


[31] 


CRESST II 


CaW04: 10.448 


10 


100 


Ua+O, 15-40 keV: <6 events 
W, 12-40 keV: <2.3 events (f) 


[32] 


DAMA/Xe-2 


Xe: 1763.2 


30 (0) 


100 


it) 


[33] 


DAMA/Nal 


Nal: 107731 


I: 22 (0) 
Na: 6.7 (0) 


100 


2-6 keVee: 0.0200±0.0032/kg-day-keVee (•) 


m 


Elegant V 


Nal: 111854 


I: 44 (0) 
Na: 13.4 (0) 


100 


4-5 keVee: 0.009±0.019/kg-day-keVee (•) 


[34] 



TABLE I: Experimental constraints used in this study. Notes to the table: (f) upper limit assuming no detected event; (J) 
limits generated based upon the binned data available in the corresponding papers (does not include any background subtraction 
and may be conservative),; (o) from an electron equivalent threshold of 13 keVee, using the quenching factor Q — Ecc/E equal 
to 0.44 for Xe [3^ : (o) from an electron equivalent thresholds of 2 keVee (DAMA/Nal) and 4 keVee (Elegant V), using Q equal 
to 0.09 for I and 0.3 for Na jsj; (•) amplitude of annual modulation. 



as follows: 7 < corresponds to an ellipse, 7 > to 
a hyperbola, and 7 = to two parallel lines. Due to 
our conies always being centered at the origin, there is a 
symmetry under (ap,a„) — > (— Op,— a„). 

Null result experiments provide an upper bound Nmax 
on the number of recoils or modulation amplitude. This 
upper bound restricts the allowed ap-a„ parameter space 
to the region inside the ellipse defined by Nmax (for the 
case in which Eqn. Hll|l describes an ellipse), or to the 
region containing the origin between the two branches 
of the hyberbola or the two parallel lines defined by 
Nmax (for the case in which Eqn. (|ll(l describes a hy- 
perbola or two straight lines). Couplings outside these 
regions would lead to more recoils than observed; cou- 
plings within these regions yield less recoils than Nmax- 

Experiments which obtain positive signals (such as 
claimed by DAMA) restrict the ap-a„ parameter space 
to a band in the shape of the appropriate conic. The 
width of the band is defined by the uncertainty in the 
observed signal (the uncertainty and hence the width of 
the band can be given to various degrees of accuracy, 
e.g. 1 or 2a). For example, if the conic is a circle, the al- 
lowed region would be an annulus; regions "outside" this 
annulus would give either too small or too large a signal. 

The relative magnitudes of A, B, and C in Eqn. 111|) 
are primarily dependent upon the structure functions of 
Eqn. I|10|l . At zero recoil energy, the structure functions 
in Eqn. (|11|) for a single nuclear element form a complete 
square, of the form N cx {aap + (3an)^ (see, e.g., Eqn. (13) 
of Ref. 123). This is a rotation of Eqn. Hence for 

g = 0, we have 7 = and the appropriate conies are two 
parallel lines. However, at nonzero momentum transfer, 
interference and spin effects destroy the symmetry so that 
different conies may result. In addition, detectors with 
multiple elements, such as the Nal of DAMA, have struc- 
ture functions from all of these elements contributing to 
the coefficients in Eqn. I|ll|l ; the resulting conic may be 



very different to that obtained from the individual ele- 
ments (see e.g. Figure 0J. For the experiments studied 
in this paper, we find that the relevant conies are either 
two parallel lines or ellipses. Some of the ellipses are so 
elongated as to appear as two parallel lines in our plots. 

We will examine the allowed shapes in the ap-a„ plane 
for current experiments and then determine where they 
overlap. In this way we will find those regions in pa- 
rameter space that are consistent with all the current 
experimental results. 



V. EXPERIMENTS 

The constraints from different experiments on the Up- 
a„ phase space are highly dependent on the choice of 
detector materials. Odd Z elements, such as Sodium, Io- 
dine, and Hydrogen, are more sensitive to Up and thus 
generate narrow ellipses with the semi-minor axis close 
to the Op axis (interference effects can cause the ellipse 
to be rotated slightly) . Conversely, elements with an odd 
number of neutrons, such as Xenon-129, Silicon-29, and 
Germanium-73, are more sensitive to and their as- 
sociated ellipses have semi-minor axes near the an axis. 
Spinless elements such as most isotopes of Silicon and 
Germanium (used by many experiments looking for spin- 
independent couplings), provide essentially no constraint 
for either ap or a„. 

Tabled (extended from (^]) displays the various exper- 
imental constraints used to generate our results. In Fig- 
ures ^fe [3 we show the limits from various experiments 
on the WIMP-proton and WIMP-neutron allowed cou- 
plings for the cases a„ = and Up = 0, respectively. Due 
to the novel detection technique used by SIMPLE (and 
the additional complexities involved) [33, we did not re- 
analyze their data in our study, but include in Figure [21 
the WIMP-neutron cross-section presented in Figure 1 
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FIG. 1: WIMP-proton cross-section limits for the case a„ = 0. The CDMS Si+Ge hmit overlaps that of Ge only; the Si is 
relatively insensitive in this case. Super-K only analyzed their data for WIMP masses above 18 GeV; their limit is taken from 
Figure 14 of Baksan analyzed fluxes down to 13 GeV. Super-K and Baksan rule out the DAM A results over their analysis 
ranges and CRESST I limits DAM A at low masses, but WIMPs between 6-13 GeV are consistent with all results for this case. 



of 0- The full ZEPLIN I results were not available dur- 
ing the writing of this paper so were not included in our 
analysis; however, a preliminary limit presented by the 
ZEPLIN group is also included in Figure 12 . 

Several experiments providing recent results, includ- 
ing Elegant V (Nal), CRESST I (Al), SIMPLE (F), 
and DAMA/Nal, used odd Z nuclei and were thus more 
sensitive to the WIMP-proton coupling ap. Other ex- 
periments, including DAMA/Xe-2 (Xe-129), CRESST II 
(Ca-43, W-183, and 0-17), Edelweiss (Ge-73), ZEPLIN I 
(Xe-129 and Xe-131), and CDMS (Si-29 and Ge-73) used 
neutron odd nuclei and were more sensitive to the WIMP- 
neutron coupling a„. We note, however, that each of 
these experiments has some sensitivity to both ap and a„, 
as the minor axes of their ellipses are not exactly aligned 
with either Up or a„. Hence, surprisingly, CRESST and 
SIMPLE obtain very good results on WIMP-neutron cou- 
pling due to their low thresholds. 

The only experiment to claim a result, DAMA/Nal, 
uses odd Z Nal detectors and a large exposure to 
search for an annual modulation. They observe 
an annual modulation amplitude of 0.0200 ± 0.0032 
events/kg/day/keVee over electron-equivalent recoil en- 
ergies of 2-6 keV. Such a modulation has been shown 
to be generally inconsistent with other experiments for 
spin-independent interactions, but we shall see there are 
choices of parameters still allowed for the spin-dependent 
case. The DAMA limits we place in ap-a„ space are 
those that reproduce the above modulation amplitude, 
using the structure functions given by Ressell and Dean 



[3^ . The WIMP cross sections on protons and neutrons 
implied by the DAMA result are plotted in Figures 
& 121 respectively. Kinematics alone can restrict the al- 
lowed masses of this modulation: small masses do not 
generate enough recoils above threshhold. One can see 
this in the figures as DAMA begins to lose sensitivity 
below about 4 GeV, requiring large cross-sections to ac- 
count for the observed signal. This feature is common 
to all dark matter experiments and lower threshold ener- 
gies are required to detect the smaller recoil energies of 
low masses. This lower limit on sensitivity is also depen- 
dent on the parameters of the velocity distribution and 
may change significantly in the presence of streams p|. 
Larger masses are also ruled out as the phase of the mod- 
ulation would be reversed. For the 2-6 keV data plotted, 
the limit is roughly 180 GeV; however, this limit is de- 
pendent upon the range of recoil energies. Freese and 
Lewis have examined this phenomenon and determined 
that masses above 103 GeV would reverse the modula- 
tion at 2 keV for DAMA In addition to the 2-6 
keV data, DAMA released modulation amplitudes for 
slightly different recoil bins, namely for recoil energies 
2-4 keV (0.0233 ± 0.0047 /kg/day/keVee) and 2-5 keV 
(0.0210 ± 0.0038 /kg/day/keVee). At the masses we are 
considering (< 100 GeV), the results of this paper are 
the same for any of the three DAMA binnings. 

All other experiments have null results which place 
bounds on the couplings. Elegant V used Nal detec- 
tors to search for an annual modulation and found no 
counts above the statistically limiting background [s^ . 
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FIG. 2: WIMP-neutron cross-section limits for the case Op — 0. SIMPLE limits are taken from Figure 1 of 0- The addition 
of the Si data clearly benefits CDMS at lower masses, ruling out DAMA for this simple case. 



While providing one of the best exclusion limits for a 
proton-only coupling at low WIMP masses, the sensi- 
tivity of this experiment is not enough to examine the 
DAMA observed region (see Figure P). CRESST I 
and SIMPLE 35] provide similar limits on WIMP-proton 
coupling. As mentioned above, due to their very low 
threshholds, CRESST I and SIMPLE also provide what 
have been the best neutron-only limits at low WIMP 
masses. DAMA/Xe-2, a neutron-odd experiment bounds 
the WIMP-neutron cross-section above 40 GeV [33l| . 
Our examination of Edelweiss (similar to CDMS, de- 
scribed below), shows that they produce a similar limit 
to DAMA/Xe-2 [s^ (see Figure llj. 

Super-Kamiokande, on the other hand, is an indirect 
detection experiment, searching for high-energy neutri- 
nos produced by the annihilation of WIMPs in the Sun's 
core after being gravitationally captured. Since the Sun 
is predominantly composed of hydrogen, the capture 
rate (and thus neutrino flux) is particularly sensitive to 
the WIMP-proton couphng Up. The Super-K detector 
searched for this additional neutrino flux by observing 
upward moving muons (induced by muon neutrinos) [14| . 
The lack of an additional muon signal leads to the most 
significant Up bounds above the analysis threshold WIMP 
mass of 18 GeV; below this mass, a significant portion of 
the muons would stop in the detector, rather than pass 
through, making the data more difficult to analyze. Dis- 
continuities of the Super-K published limits occur at 80 
GeV (W mass) and 174 GeV (top mass) due to differ- 
ent possible annihilation products. The Super-K WIMP- 
proton cross-section limits are shown in Figure ^ these 



limits are taken directly from Ref. |l4j . derived therein 
using the model-independent technique of Kamionkowski 
et al. [S^. We note that the bounds from Super-K rely 
on three assumptions: First, we assume that the WIMPs 
have achieved equilibrium in the Sun, so that the capture 
rate is equal to the annihilation rate; this assumption 
is almost certainly correct, particularly at small WIMP 
masses. Second, we assume that either WIMPs are equal 
to their antipartners so that they can annihilate among 
themselves or, if WIMPs are not the same as the anti- 
WIMPs, that there is no WIMP/antiWIMP asymmetry. 
Third, the WIMPs do not decay to light fermions only 
(in which case fewer observable neutrinos would be pro- 
duced) . The second and third ways to evade the SuperK 
bounds were pointed out by Kurylov and Kamionkowski 
jS]. Here we assume that all three assumptions are cor- 
rect; then Super-K rules out the proton only coupling for 
WIMP masses above 18 GeV. 

Baksan, another indirect detection experiment, pro- 
vides the tightest bounds on proton only coupling be- 
tween 13 and 18 GeV. Although Baksan provides weaker 
WIMP-induced muon flux limits than Super-K above 18 
GeV, the Baksan experiment has analyzed these fluxes 
down to a lower WIMP mass of 13 GeV [l^. These 
flux limits can be used to determine WIMP-proton cross- 
section limits in the same manner as Super-K; such an 
analysis is not available from the Baksan collaboration at 
this time. However, the Baksan cross-section limit can be 
determined by rescaling the Super-K limit by the relative 
flux limits of these two experiments (since the expected 
flux is proportional to this cross-section). Below 18 GeV, 
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where there is no Super-K hmit to rescale, we extrapo- 
late the Baksan cross-section hmit using their given flux 
hmit and the WIMP mass dependence of the expected 
flux, given by Eqn. (9.55) of Ref. "2^. We find that Bak- 
san rules out proton only coupling down to WIMP masses 
of 13 GeV. 

CDMS II, in the Soudan Underground Laboratory, 
uses Silicon and Germanium detectors with strong dis- 
crimination to provide the strongest spin-independent 
cross-section limits to date |23|- Analysis of CDMS re- 
sults in the spin-dependent sector have consistently been 
ignored, as the small natural abundance of spin odd iso- 
topes in Sihcon (4.68% Si-29) and Germanium (7.73% 
Ge-73) has given the impression that CDMS is not sen- 
sitive enough for any significant spin-dependent results. 
Our analysis, however, shows that the extremely clean, 
background-free data is sufficient to overcome the smaller 
spin sensitive mass, with the odd neutron Ge-73 (Silicon 
data has not yet been published) providing significant 
limits to Qn- The limits we obtain are due to the fact 
that no events were observed in 52.6 days of livetime for 
four 0.250 kg Germanium detectors over 10 keV to 100 
keV recoil energy, using an efficiency of 0.228 for 10-20 
keV and 0.38 for 20-100 keV; an efficiency of 0.076 for 
5-10 keV is used to determine limits at lower thresholds. 
Since Silicon is lighter than Germanium, it will provide 
better limits at small WIMP masses. As the Silicon data 
has not yet been published for the first run of CDMS II, 
we augment the CDMS II Ge results with the most recent 
CDMS I Si data The last CDMS I run used the same 
detector tower as now used in CDMS II, containing two 
0.100 kg Silicon detectors as well as the four Germanium 
ones. However, one of these Si detectors was contami- 
nated and not used in the data analysis. No Si events 
were found in the range 5-55 keV after 65.8 days of live- 
time, with the same efficiencies as above. While several 
Si events were observed above 55 keV, these are consis- 
tent with a neutron background. Even if these events 
were WIMP scatters, they are not consistent with a light 
mass. Since the Si data is being used to extend the limits 
at low masses (Ge data dominates at higher masses), we 
feel justified in ignoring the high energy events. CDMS 
hopes to install a total of five towers and achieve a total 
exposure during the CDMS II run of 1200 kg days for 
Ge and 300 kg days for Si; we will use these numbers to 
project future limits, assuming a 5 keV threshold and null 
results [2^. We have used the structure functions given 
by Ressell et al. for Si-29 and Dimitrov et al. for 
Ge-73; a comparison of the Ge-73 functions from Ressell 
shows the results are relatively insensitive to the models 
used to derive the form factors. 



VI. RESULTS 

We first present results for the simple cases where the 
WIMP couples to only the proton or neutron. From a 
combination of constraints in Figure [21 one can see that 



DAMA results are not compatible with limits from other 
experiments for the case of the WIMP coupling with 
neutrons only. However, one can see in Figure HI that 
DAMA still survives for the case of proton only coupling 
for masses less than 18 GeV. 

Figure ^ displays the current limits for the WIMP- 
proton cross-section assuming the WIMP couples only to 
the proton (a„ = 0). CDMS does not currently constrain 
this coupling. The Super-K and Baksan results, however, 
strongly rule out the WIMP-proton coupling that would 
be required to explain DAMA for masses above 13 GeV. 
Elegant V still allows for the DAMA observed modula- 
tion at these lower masses and CRESST I provides only 
a lower mass limit of 6 GeV. 

Figure [21 displays the WIMP-ncutron cross-section in 
the alternative case that the WIMP couples only to the 
neutron Op = 0. DAMA/Xe-2 and Edelweiss rule out 
this case for masses above 20 GeV. The surprising result 
is that CDMS II significantly improves upon this limit, 
ruling out masses above 10 GeV. Furthermore, if one in- 
cludes the CDMS I Si data, CDMS rules out aU of the 
DAMA allowed region for this case. Recent ZEPLIN re- 
sults are extremely powerful as well. 

While the above two cases demonstrate the different 
sensitivities of the experiments, they are not necessar- 
ily representative of the more general parameter space, 
where both ap and a„ may be non-zero. Figure [3 shows 
the Up-Qn regions allowed by DAMA, CDMS, and a com- 
bination of the remaining experiments at several different 
masses for this more general case. 

The DAMA observed modulation generates an ellipti- 
cal ring, flattened in roughly the direction of the Op-axis 
(due to the odd Z Nal). [At low masses, the portion of 
the elliptical ring that is plotted looks like two parallel 
bands.] The ellipse is slightly rotated from the ap-a„ 
axes; this rotation increases at larger masses as observed 
recoils tend to come from the Iodine atoms rather than 
the lighter Sodium (which has different structure func- 
tions). The DAMA allowed couplings are those that fall 
on the elliptical bands, not between them; couplings "in- 
side" the inner edge of the elliptical ring (between the 
two bands of the ellipses) are too small to generate the 
observed modulation. 

Other experiments, which have null results, place 
bounds on the ap-a„ parameter space, as shown in Fig- 
ure [3 The light grey (pink) horizontal regions, labelled 
as "other" in the legend of the figure, illustrate a com- 
bination of limits from Super-K, Baksan, Elegant V, 
CRESST, Edelweiss, and DAMA/Xe-2 (we remind the 
reader that DAMA/Xe-2 is entirely different from the 
DAMA/Nal annual modulation results). Only couplings 
within these regions satisfy the null result constraints of 
all these experiments. Super-K and Baksan results de- 
pend upon WIMPs interacting with hydrogen, and thus 
inherently limit the proton coupling ap only. The Super- 
K and Baksan allowed regions are bands along the a„- 
axis (corresponding to the region between two parallel 
lines). The DAMA/Xe-2 and Edelweiss allowed ellipses. 
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FIG. 3: Spin-dependent allowed couplings at several different WIMP masses (in different panels). The legends in the figure 
indicate regions allowed by different experiments. The thin black horizontal bands (ellipses at higher masses) show the region 
allowed by DAMA (note that the black bands are indicated as white for regions allowed by all experiments); the allowed 
couplings for DAMA are those that fall on the bands, not between them. The dark (blue) vertical regions correspond to 
couplings allowed by the CDMS I Si and CDMS II Ge data sets. The remaining (pink) horizontal region, labelled as "other" in 
the legend, represents couplings allowed by the combined null results of Super-K, Baksan, Elegant V, CRESST, DAMA/Xe-2, 
and Edelweiss (couplings outside of this region are ruled out by at least one of these experiments). The region allowed by all 
experiments is indicated by white bands. Only inside these white bands is the DAMA signal consistent with the null results 
from other experiments. For masses above 13 GeV there is no consistent region (no white band). 



whose minor axes are along the a„ direction, are far 
broader in the Op direction. The intersection of these 
regions (horizontal Super-K & Baksan bands with ver- 
tical DAMA/Xe-2 & Edelweiss ellipses) forms a roughly 
rectangular region above 13 GeV, clearly visible in the 
15 GeV panel and getting smaller at higher masses. Any 
allowed model must lie within this "rectangle." Below 
18 GeV, Super-K and DAMA/Xe-2 no longer provide 
the best limits (or any limit, in Super-K's case); below 
13 GeV, Baksan and Edelweiss also lose their sensitivity. 
For these lower masses. Elegant V and CRESST I (in- 
cluded in the region labelled "other") provide the best 
limits; however, they do not significantly constrain the 
DAMA space between 5 and 13 GeV. 

The recent CDMS results are also shown in Figure 13 
CDMS's null result, based upon odd neutron Si and Ge, 
generates an allowed region elongated along the a^-axis. 
In FigureEl the region allowed by the CDMS data is dark 
grey (blue) and is labelled "CDMS I&II." The addition 
of CDMS I Si data increases sensitivity at low masses, 
due both to the lighter mass and lower threshhold than 
the CDMS II Ge data. As the sensitivity changes from 



predominantly Ge at high masses to Si at low masses, the 
shape of the allowed region rotates slightly (as Si and Ge 
have different structure functions). 

In sum, we have applied null results from all other 
experiments to see whether or not DAMA annual modu- 
lation results can be compatible. The only remaining re- 
gions are shown in white bands in FigureEl Only WIMP 
masses in the range 5-13 GeV survive as spin-dependent 
candidates matching all the data. In the next section 
(Summary and Discussion) we discuss the viable ranges 
of couplings in the ap-a„ plane. 

Note that, when one allows both ap and a„ to be 
nonzero, the bounds of Figures ^ & [3 on proton-only 
and neutron-only couplings can be considerably weak- 
ened. We illustrate this with an example in Figure 01 
which shows the limits in the ap-a„ plane separately for 
the CDMS I Si and CDMS II Ge at a WIMP mass of 10 
GeV. While Ge limits Op to below 160 for a„ = 0, even 
a small non-zero value of a„ allows for ap values several 
times larger (e.g. we can have ap = 330 at a„ = —20). 
The CDMS II Ge limit in FigureQlfor the WIMP-proton 
coupling applies only to the case of a„ = 0, but the limit 
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10 GeV 



CDMS I Si 
CDMS II Ge 
Combined 



FIG. 4: Spin-dependent couplings allowed by CDMS at a 
WIMP mass of 10 GeV. The separate limits from Si and Ge 
each produce extremely long ellipses; the combined limit is 
significantly smaller. 



in the general case is more than an order of magnitude 
weaker. 

Referring again to Figure^] the Si data forms a similar 
elliptical limit as the Ge, but slightly rotated. Even the 
slight rotation is significant enough that the two narrow 
ellipses rarely overlap- the combined Si and Ge limit is 
significantly smaller than that of either of the separate 
limits. Detectors employing multiple elements can break 
the (near) degeneracy typical of individual elements. 



VII. 



SUMMARY AND DISCUSSION 



Figure 131 illustrates our main results. When we allow 
both Up and a„ to be nonzero, there is no region of pa- 
rameter space with WIMP masses above 13 GeV which 
is compatible with both the positive results of DAMA 
and the null results of Super-K, DAMA/Xe-2, Edelweiss, 
and CDMS. Below 13 GeV, there are line segments in the 
ttp-cLn plane which are still compatible with the DAMA 
data and all null results from other experiments; CDMS 
II is the most constraining. At 12 GeV, the compat- 
ible line segments are ap — — 0.084a„ -I- 1.94 (±0.16) 
with -4.62 < a„ < 4.39; at 8 GeV they are Up = 
-0.084a„ + 2.81 (±0.23) over -15.8 < a„ < 15.8 (these 
describe only one of the line segments, the other is found 
by taking (a„,ap) (— a„,— Op)). This includes cases 
for which the proton and neutron couplings are compa- 




rable (e.g. 



= a„ = 1.8 at 12 GeV and 



2.6 



at 8 GeV), in which case the WIMP-proton and WIMP- 



FIG. 5: Spin-dependent allowed couplings at 5 GeV for 
DAMA (black elliptical band), CDMS (vertical region), and 
CRESST I (angled region). The CDMS and CRESST I 
null results only allow for couplings within the intersection 
of their corresponding allowed regions (hatched region); this 
intersection does not include any couplings that can repro- 
duce the DAMA observed modulation. Below 5 GeV, the 
CDMS/CRESST I constraint grows rapidly stronger relative 
to the DAMA required couplings. Note the axes are not at 
the same scale. 



neutron cross-sections are of the same order (equal for 
the given examples). We also note that only models with 
jflpl > 1 are consistent at any mass. Below about 5 
GeV, CRESST I becomes significant enough to exclude 
(in combination with CDMS) all of DAMA's parameter 
space (Figure |SJl. 

We note that the small remaining regime of spin- 
dependent parameter space could be confirmed or ruled 
out by more complete analysis of Super-K or other so- 
lar neutrino (indirect detection) data. If Super-K, Bak- 
san, and ANTARES (or in the future IceCube and 
ANTARES) were able to push down their analyses be- 
low 13 GeV WIMP masses, they would be able to ei- 
ther find the WIMP annihilation from the Sun compat- 
ible with DAMA, or rule out DAMA entirely as due to 
spin-dependent interactions of a WIMP that is its own 
anti-WIMP. If the WIMP and anti-WIMP abundances 
in the Sun would differ substantially, there would be no 
annihilation signal from the Sun even if WIMPs may be 
abundant enough to give rise to a signal in DAMA. Pre- 
viously UUio et al. 5] attempted to push the Super-K 
results down to lower WIMP masses. However, we feel 
that their analysis may have been too constraining due to 
the fact that the energy threshold for the muons in the ex- 
periments was not taken into account (it was set to zero). 
We feel that this analysis must be left to the experimen- 
talists. We encourage the experimentalists to study the 
bounds at lower WIMP masses that would arise when 
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one extends the analysis to more difficult cases: muons 
that stop in the detector, muons that start in the de- 
tector, and contained events (muons that start and stop 
in the detector). For the contained events one should 
consider both muon and electron neutrinos. If such an 
analysis were possible, then one could confirm or refute 
the remaining parameter region for spin-dependent inter- 
actions in DAM A. 

In the future, as other experiments become more sen- 
sitive, they will of course further restrict the parameter 
space as well. New CRESST I or SIMPLE data should 
also be able to push down the WIMP mass compatible 
with DAMA. Edelweiss, CRESST II and ZEPLIN [13, 
whose role is similar to CDMS in this regard (their el- 
lipses align with the ap axis), will also become important. 

Accelerator bounds on spin-dependent interactions are 
far weaker than the limits we have discussed from dark 
matter experiments. At the high energies in accelera- 
tors, the WIMPs couple directly to the quarks and can 
no longer be treated as effectively coupling to the nu- 
cleons as a whole. Such scattering may occur through 
exchange of several different particles, leading to interfer- 
ence effects that further reduce the cross-section. Likely 
signatures of such scatters are also either unobservable or 
dwarfed by backgrounds. Looking at the case of a simple 
Z exchange, the lowest order case oi qq ^ Z + X with Z 
producing a WIMP pair Z — > XX would not present an 
observable signal as the WIMP itself, being electrically 
neutral and weakly interacting, is not directly detectable. 
The signal is potentially observable if a gluon is also emit- 
ted, qq ^ Z + g + X, but other Standard Model processes 
produce the same signatures at rates much larger than 
those expected for Up and a„ of order one. Even assum- 
ing no interference effects, the experimental limits from 
this and similar processes give Op < O(IO^), much weaker 
than the limits of 0(1) from the dark matter experiments. 
As a final remark, we note that the the Minimal Su- 
persymmetric extension of the Standard Model does not 
provide a dark matter candidate with the characteristics 
pointed to by our study. 
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the nucleus recoil energy to find the recoil rate R per unit 
detector mass: 



EilQ 
Ei/Q 



dEe{E)-^a{q)7^{E,t). (A.la) 



e{E) is the (energy dependent) efficiency of the exper- 
iment, due, e.g., to data cuts designed to reduce back- 
grounds. Q is the quenching factor relating the observed 
energy Edet (in some cases referred to as the electron- 
equivalent energy) with the actual recoil energy E^-ec- 
Edet = QErec- The energy range between Ei and E2 
is that of observed energies for some data bin of the de- 
tector (where experiments often bin observed recoils by 
energy). Note some single element detectors can cali- 
brate their energy scales to Edet = Erec, in which case 
the quenching factor can be ignored. For detectors with 
multiple elements, the total rate is given by: 



i?tot(t) = ^/^i?^(<) 



(A.lb) 



where fi is the mass fraction and Ri is the rate 
(Eqn. ljA.la|) ') for element i. 

The expected number of recoils observed by a detector 
is given by: 



Nrec = MdetTR 



(A.2) 



where Mdet is the detector mass and T is the exposure 
time. To calculate the constants in Eqn. Hll() . we first 
rewrite Eqn. H10|) as: 

o-soiq) = hpp{q)al + hpn{q)apan + /i„„(g)a^ (A. 3) 

where 



h^j{q) 



2J+1 



■Stj{q) 



(A.4) 



is independent of the couplings Up and a„ . The constants 
Arec, Brec, ^^d Crec in Eqn. (|lla|l are then determined 
from Eqn. IIA.2p by replacing a{q) in Eqn. ljA.la|) with 
the appropriate hij{q), e.g. 
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APPENDIX: DETECTOR RECOILS 



E2/Q 
Ei/Q 



Arer = MdetT / dE e{E) hpp{q) 7^{E, t). 



(A.5) 

Similarly, B^ec is given by Eq. IjA.Sp but with hpp replaced 
by hpm and Crec is given by Ea. (|A.5p but with hpp re- 
placed by hnn- 

For large detectors that can obtain high statistics, the 
annual modulation of the rate (due to the rotation of the 
Earth about the Sun) may be detectable, where 



R{t) w Ro + R,nCos{ujt). 



(A.6) 



To determine the number of expected recoils for a given 
experiment and WIMP mass, we integrate Eqn. Q over 



The rate (for most energies and for a Maxwellian velocity 
distribution) is maximized around June 2 and minimized 



11 



around December 2; this phase is reversed at high enough 
energies |33- The modulation ampUtude, 

Nma = Rra^\ [i?(Junc 2) - i?(Dcc 2)] , (A.7) 

denoted by Nma for consistency with Eqn. Ullafl . can be 
put in the form of Eqn. (|llb|l in the same manner as the 
total recoils above. Hence A^aa is the value of Eqn. ljA.7p 
when a{q) hpp{q) in Eqn. (jA.lajl : 

rE2/Q 

Ama = / dEe{E)-^hpp{q) 

JEi/Q ^"^M 

X i[?7(£;, June 2) - 77(£;,Dec 2)].(A.8) 



Similarly, B^a is given by Eq. ljA.8|) but with hpp replaced 
by hpn, and Cma is given by Ea. l)A.8|l but with hpp re- 
placed by hnn- 
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